Tau functions for the Dirac operator on the cylinder 



O. Lisovyy *> t 



* Bogolyubov Institute for Theoretical Physics 
Metrolohichna str., 14-b, Kyiv-143, 03143, Ukraine 

t Departement de Mathematiques, Universite d'Angers, 
2 Boulevard Lavoisier, 49045 Angers, France 

Abstract 

The goal of the present paper is to calculate the determinant of the Dirac operator with a 
mass in the cylindrical geometry. The domain of this operator consists of functions that re- 
alize a unitary one-dimensional representation of the fundamental group of the cylinder with 
n marked points. The determinant represents a version of the isomonodromic r-function, 
itroduced by M. Sato, T. Miwa and M. Jimbo. It is calculated by comparison of two sec- 

■ tions of the det*-bundle over an infinite-dimensional grassmannian. The latter is composed 
of the spaces of boundary values of some local solutions to Dirac equation. The principal 

■ ingredients of the computation are the formulae for the Green function of the singular Dirac 
operator and for the so-called canonical basis of global solutions on the 1-punctured cylinder. 
We also derive a set of deformation equations satisfied by the expansion coefficients of the 
canonical basis in the general case and find a more explicit expression for the r-function in 
the simplest case n = 2. 
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1 Introduction 

The main objective of quantum field theory is the calculation of correlation functions of local 
operators, usually represented via functional integrals 



For a generic interacting QFT such calculation can be done only by means of perturbation theory. 
^ ■ However, in two dimensions there is an interesting way to construct an interacting theory from 

the free one. Let us consider the action of free massive Dirac fermions in the flat spacetime, 



S[il>,ij)] = -j d 2 x^D^. 
Correlation functions of the (interacting) monodromy fields are defined as 

The integration in the numerator is performed over field configurations that are branched at 
the points a v with the monodromies e 2mXjy (u = 1, . . . , n). This change of domain of the Dirac 
operator is symbolically reflected by indexing D a,x . The integrals in (II. ID can be formally 
evaluated to the determinants of the corresponding operators, 

det D a,A 

(O x Hai)...0^(a n )) = ^-. (1.2) 
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Note, however, that the RHSs of both (jl.lj) and (|1.2|) are equally ill-defined quantities. 

The determinants of Dirac operators on compact manifolds are usually determined via the 
(^-function regularization. Starting from [16], they have been extensively studied in the mathe- 
matical literature. The massless Dirac operators on Riemann surfaces deserved special attention, 
as multiloop contributions to the partition function in the string theory are expressed through 
their determinants (rigorously defined by D. Quillen in [15]). These determinants can be thought 
of as the functions on the moduli space of complex structures on the surface. 

In the case we are interested in the Dirac operator is defined not on a compact manifold, but 
on the universal covering of a surface with marked points, the determinant being a function of 
their positions. The problem of rigorous definition of the determinant and the Green function 
for the Dirac operator with branching points on the Euclidean plane was solved by Palmer in 
|10j . His work relies heavily on the analysis of monodromy preserving deformations for the 
Dirac operator, developed earlier by Sato, Miwa and Jimbo [17]. More precisely, Palmer's 
determinant represents another version of SMJ r-function. Its logarithmic derivatives with 
respect to the coordinates of branching points are expressed via the expansion coefficients of 
some special solutions to Dirac equation, that can be constructed from the so-called canonical 
basis of solutions. The theory of isomonodromic deformations gives a set of nonlinear differential 
equations satisfied by these expansion coefficients. Moreover, in the simplest case n = 2 an 
explicit formula for the determinant was found |10j . 

Later similar results were obtained for the Dirac operator on the Poincare disk [8j [12j [EH]. 
In this connection we should also mention the recent work of Doyon [5], where the two-point 
correlation function of monodromy fields in the hyperbolic geometry was calculated by field- 
theoretic methods. 

In the present paper, we define and calculate the determinant of the massive Dirac operator 
in cylindrical geometry. The latter corresponds to QFT in the finite volume or at non-zero 
temperature. This work was inspired by recent progress in the study of Ising model — calculation 
of finite-size correlation functions [U [2], spin matrix elements [3] and direct derivation of the 
differential equations satisfied by two-point correlator in the continuum limit [7J . Ising model is 
related to the above theory with special monodromy \ v = 

This paper is organized as follows. In the next section we introduce the canonical basis 
of global solutions to the Dirac equation on the cylinder and calculate it explicitly for n = 1 
(see Theorem I2.3f) . All subsequent computations are based on these formulae. The Green 
function of the singular Dirac operator is defined in Section 3. Its derivatives with respect to 
the coordinates of branching points are expressed through some solutions to the Dirac equation 
and have remarkable factorized form (formulae (|3.23[) and (|3.24p ). At the end of the section, 
the Green function on the 1-punctured cylinder is computed (see (I3.29p - (l3.30j) ). Section 4 is 
devoted to the definition and calculation of the r-function. We introduce the det*-bundle over 
an infinite-dimensional grassmannian that consists of the spaces of boundary values of some 
local solutions to the Dirac equation. The r-function is obtained by comparison of the canonical 
section of this det*-bundle with a section that is constructed using the one-point Green functions. 
The logarithmic derivatives of the r-function can also be written in terms of the expansion 
coefficients of the solutions (j3.20p . This shows that it is independent of the chosen localization. 
As an illustration, we find a more explicit expression for the r-function when n = 2 (formula 
(|4.25p ). Finally, in the Section 5 a set of deformation equations for the expansion coefficients 
is derived. We conclude with a brief discussion of possible generalizations, open problems and 
application of obtained results in the quantum field theory at non-zero temperature. 
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2 Canonical basis of solutions to Dirac equation 



2.1 Definitions 

Let a = (ai, . . . ,a n ) be a collection of n distinct points on the cylinder C. The fundamental 
group 7Ti (C\a; xq) is generated by homotopy classes of n + 1 loops 7o, • • • > 7n shown in the Fig. 1. 
It acts on the universal covering C\a by deck transformations. Let us fix a one-dimensional 
unitary representation 

p x : TT 1 (C\a;x )^U(l), [ 7l/ ] ^ e - 27riA ", i/ = 0,...,n, (2.1) 

A gM, A„eR\Z, i/ = l,...,n. 




Fig. 1. 



As usual, we replace the cylinder by the strip S = {(x, y) G M 2 : < y < (3} whose upper 
and lower edges are identified. The Dirac operator on C\a is induced by the Dirac operator on 
M 2 , which can be written as 

D = (-l "|")' (2 ' 2> 

where z, z — standard complex coordinates 

z = x + iy, \d z = ~(d x -id y ), 

z = x-iy, [dz = \ (d x + id y ). 

We are looking for multivalued solutions tp : C\a — ► C 2 to Dirac equation that transform 
according to the representation (|2.ip . 

Dip(x)=Q, t/j^x) = p\{[i\) ■ ip(x). 

This problem can be reformulated as follows. Fix a system of branchcuts b — • • • 3 ^nj 
do, . . . , d n ) shown in the Fig. 2 and consider the solutions to Dirac equation on C\b that can be 
continued across the branchcuts away from the points a±, . . . , a n . The solutions we are interested 
in have left and right continuations across b v that differ by the factor e 2nlXv (y = 1, . . . , n). The 

V 

continuations across d v differ by exp(27ri ^ A&) , v = 0, . . . , n. 

k=0 
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Fig. 2. 



To describe the local behaviour of such solutions in the neighbourhood of the point a u , 
consider an open disk B of sufficiently small but finite radius, centered at a v , and introduce in 
B polar coordinates 

1 1/2 ( f) — I e -iH> 



zi z—a v ' 
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The local form of the Dirac operator on B is then 

1 / 171 —( 

2 



D 



-up 



(Or 
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Since for any multivalued solution tp the function e~ lXiy,P/ ip is single-valued on B, it can be 
expanded in Fourier series. Substituting the series into Dirac equation, one obtains [17] 

Yl { a k w k+\Ao-u\ +b k wZ_ Xv [a v \} . (2.3) 



where 



wi [a u \ 



e iG-l/2)¥'J z _ 1/2 ( mr ) 
e i( - l+1 / 2 ^Ii^/9.{mr) 



w l [a„ 



(2.4) 



e -i(i+-imv>i l+1/2 ( mr ) 

H+1/2V"') J V e-^-^I^imr) 

and Ii(x) is the modified Bessel function of the first kind. 

To obtain some kind of regularity, we put certain conditions on the singular behaviour of the 
function ip at the point a u . There are two essential types of constraints: 

• Let < Ajy < 1 and require t/j to be square integrable in the neighborhood of a u . When 
\z\ — > a u , the asymptotics of special solutions has the form 



wi [a v \ 



( (m(z-a„)/2) 1 I \ 

F|)i ~ 

(m{z-a v )/2) l+ 1 



+ ..., 



( (m(z-a„)/2) l+ i 
(m(z-a„)/2)'-? 



where factorials are understood as V. = T(l + 1). Then to satisfy the condition of square 
integrability, a part of coefficients in (|2.3p must vanish, 



ip[a v ] = o_i/ 2 w_i/ 2+ A v [ai/] + {akW k+ \A a A + hw k *_ Xv [a v ]} . 

k>0 



(2.5) 
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• Let — \ < \ u < g and require 

where -ff 1 [a !/ ] denotes the space of functions that are single- valued and square integrable 
in the neighborhood of a v together with their first derivatives. Then 

= Yl { a kw k +\ v [du] + hwk- Xv M } • (2-7) 

k>0 

Now let us consider multivalued solutions with monodromy (|2.ip that are square integrable 
at |x| — ► oo and satisfy (12, 5p or (12. 7h in the neighborhood of each singularity. The spaces of 
solutions of the first and second type will be denoted by W a ' A and W a ' A respectively. 

Theorem 2.1 dim W a,A < n; dimW a ' A = 0. 

■ Consider a positive definite scalar product on W a,A : 

if i 
m / . . m 



(u, w) = (w, u) = J u ■ w idz A dz = J (u\Wi + U2W2) idz A dz. (2-8) 

C\a C\a 

Note that the expression under the integral is indeed a single-valued function on C\a. This 
function is integrable due to imposed boundary conditions. From the Dirac equation on C\b it 
follows that 

(s± Wl = d z w 2 , { % ux = dzU 2 , 



and we get 

Denote by D e (a u ) the disk of radius e about a v . Using (|2.9p and Stokes theorem, one obtains 

n „ 
(u,w) =im^ lim ® 
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— (u\W\ + U2W2) dz Adz = —d(u2W\dz) = d(u\W2dz). (2A 



U2W\dz 



U ~ l 8D e {a v ) 

• v I ( M 1 \ (mz/2) x » , M , , (mz/2)- x » 

im^ hm <j> la ( _l, n (u)± — J-J— + . . . + b) /Ju) —, ' , '— +■ 



'" 1 8D e {a v ) 



(mz/2) A -- 1 | (mg/2) 



1-A„ 



+ . . . ]dz 



-4 ^ 2 (u)o3 /2 H sin ttX u . (2.10) 



i/=i 

Or, analogously 



n „ 

(u,w) = -im^ lim ® u\W2dz 

U=1 dD E {a„) 
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■ v I ( H / \ (mz/2)^- 1 r v) . , (mz/2) 1 ^ 

— Ei-o / (A, -1)1 + - + b i>) (il^, + 

n 



= -4 J^ari/aHbi/iHsiiiTrA,, = («;,«). (2.11) 

If the dimension of W a ' A were greater than n, we would be able to construct a solution v G W a,A 
with all a_y 2 (v) = [y = l,...,n). This solution has zero norm (v,v) = and, therefore, 
vanishes identically, implying the first statement of the theorem. 

Note that the solutions of the second type are square integrable with respect to the inner 
product (|2.8j) . We can show in absolutely analogous fashion that (v,v) = for all v G W a,A . 
Consequently, dimW a,A = 0. ■ 

Supposqj that dim W a ' A = n. Then we can fix a canonical basis {wj^i,,,, of this space, 
having chosen a_-jy 2 (w^) = 5^: 

W M M = V!i)-l/ 2+ A„W + Yl { ai k( w n) w k+\A a v] + &( fe ) ( W M) U; jfc-Aj^]} ( 2 - 12 ) 

fc>0 

Remark. Let us calculate the inner product of two elements of the canonical basis in two ways 
- by the formula (pTTOj) and its "conjugate" (pTT]) : 



(w M ,wv) = -46^2 (w M ) sin 7rA„ = -46^y 2 (w„) sin7rA M . ( 2 -13) 

We have obtained a set of algebraic relations between the expansion coefficients by^w^. In 
what follows, we will deduce additional relations and use them in the construction of deformation 
equations. 

The "planar" analog of the previous theorem has an instructive illustration when n = 1. In 
the case of a single branching point one can suppose it to lie at zero. Then any solution with 
required singular behaviour is represented by the expansion 

if) = a_ 1/2 w_ 1/2+x [0] + {a k w k+x [0} + b k w k *_ x [0}} 

k>0 

on the whole punctured plane M 2 \{0}. This expansion will be square integrable at infinity if 
and only if 

a k = for k > 0, 
< bk = for k > 1, 

A/2 = ~ a -l/2, 

since the only integrable combinations of partial solutions (12.4j) are 

w l [0}=wl l [0}-w l [0}. 



1 The proof is based on some technique from functional analysis and is very close to the proof of Theorem 3.2.4 
in [IT]. 
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Then, as one could expect, for n = 1 the space W 0,A is generated by the single element of 
canonical basis 

W = ^_i/ 2 +a[0] - U>1/2-a[°] = -^l/2-A[0]- 

With some efforts, it is also possible to find an explicit formula for the canonical basis on the 
1-punctured cylinder. This problem will be solved in the next subsection, using a generalization 
of the method of Fonseca and Zamolodchikov [6]. 

2.2 Canonical basis on the cylinder with one branching point 

We are looking for the solution ip to Dirac equation on the strip < y < (3 

which has the following properties: 

• The continuations of this solution to the left and right halfplane are quasiperiodic in y, 

ip(x, y + p) = e 2niX °ip{x, y) for x < 0, (2.14) 
t/j(x, y + P) = e 2mX i>{x, y) for x > 0, (2.15) 

where A = Ao + Ai. 

• It satisfies the normalization condition 

lim (mz^f-^^ix^y) = , (2.16) 

i*M r(Ai) 

where the fractional power of z is defined as 

Z l-Ai =e (l-A 1 )lnz j < Im(lnz) < 2vr. 

Theorem 12.11 shows that these requirements determine the solution uniquely. 

The function e~ 2mX ° l ''^ 't/j is periodic in the left halfplane and therefore can be expanded 
there in Fourier series. Substituting the series into Dirac equation, one obtains a general form 
of the solution for x < 0, 

^<0(x,y) = -A V ^ {Bn l a e rnxcoshe n+ imysu&O n ( «M ^ 

t—* mp cosh 9 n V 1 I 

where sinh# n = and the factor 1 /(m(3 cosh 9 n ) is introduced for further convenience. Anal- 
ogously, the general form of solution in the right halfplane is 

J/) =A V f ^ e -mz^e n -imyM n ( \ (21g) 

' mp cosh # n VI/ 

neZ-A 

Of course, in order for the series (|2.17[) and (|2.18p to converge the functions G{9) and H (9) 
have not to grow too rapidly as 9 —* ±oo. Moreover, we shall assume that G(9) and H(9) are 
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> Re0 



Fig. 3. 



analytic in the strip — ^ — 5 < Im9 < ^ + 5 for some 5 > 0, so that (|2.17p and (|2,18p can be 
represented via contour integrals (see Fig. 3) 

^x<0(x,y) = A [ — . ^ 6 }. „ ., g^coshe+imysinhe ( e \ 

Y <UV ,yJ J c _ \j c+ 2lT 1 - e »m/3sinh0-27riA o V 1 / 

ih x>0 (x,y) = A f — _ e -mxcoshe-im^sinhe / ~e e \ 

Y >UV _/ c _ y c+ 27T 1 _ e -im(3sinhe~2ni\ \ 1 J 

If < y < /3, the contours C+ and C_ can be continuously deformed into ImO = ^ and 
lm9 = —5, respectively, defining the continuations of ip x< o(x,y) and i^ x >o(x,y) on the whole 



\- A f de { G(6 + iTT/2)e imxs ' mhe ~ m y cosh( 

' / 27T 1 ]^ _ g— m/3cosh0— 27riAo 



ze 17 



^<D<o(aJ,y) = ^ / — <, 5 t^oz^i; ( 1 1 + 

_ ^ 7r ^2)e _im:csinhe+mi/coshe / — ie e 

2 gm/3 cosh 0— 27riAo I ^ 



^>o(x,y) = Aj -j ^^^^ ^ i J + 

— oo 

jj^q iff / 2~) e^ mx s ' n ' 1 m ^ cosh 



^ g— m/3 cosh 9— 2wi\ V 1 

These continuations coincide if two functional relations for G{9) and -£f(#) hold 

G{9 + in/2) 1 - e~ m P cosh " 2 ™ A o 

H{9-m/2)~ 1 _ e -m/3cosh6»-27ra ' 



(2.19) 



G{e-m/2) = _ e!htiAl i_ e -^c°shQ+2^A 

H{9 + m/2) 1 _ e -m/3cosh0+27ri.A ' ^ ' ' 

The relevant solutions of these equations can be found using the following lemma. 
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Lemma 2.2 Consider two functions, f(9) and g(9), that are analytic in the strip \lm9\ < 5. If 
in this strip \f(9)\ = O f rj^rgp J and \g(9)\ = 0(1) as He 9 — > ±oo, i/ien i/ie functions 



oo oo 

/J/V /■ J/3/ 

— tanh(0' - 9)f(9'), V (9) = J — sech (9' - 9)g(9'), 9 e 



(2.21) 



can be analytically continued to the strip \ImO\ < ^ + 5. Furthermore, if |Im0| < S, these 
continuations satisfy the relations 



is 9 + 



ITT 



~2~ 



~if(0), 



ITT 



+ rj 



ITT 



9(9). (2.22) 



■ Obviously, the expressions (|2.2ip for v(9) and rj(6) are analytic functions in the strip |Im0| < 
5. Their analytic continuations to |Im0| < S + 5 are 



"(0) 



T^/(0T?)+P / ^coth(0'-0±^)/(0' 



ime=±4 



dff 



ITT 



27T 



77(0) 



^g{9^)±iPj ^-csch(0'-0±^W0' 



Imfl=±f 2 



f <Im0<f +5 



2tt 



7 7T 



?.7T 



d0' 



*7 f-^T + / — tanh(0'-0)/(0'), 



2tt 



— oo 
oo 



f <Im0<f +5 



/7T 



d0' 



5(0- y)+ / — sech(9'-9)g(9'), 



v{9) 



-%-5<lm6<-Z J \ 2 1 



-oo 
oo 



n(9) 



§-<S<Im0<-f 



Z7T 

//(0 + y) + 



2tt 
2tt 

d0' 

2^ 



tanh(0' - 0)/(0'), 
sech(0' - 9)g{6'), 



where the first two integrals are understood in the principal value sense. Then the statement of 
the lemma follows immediately. ■ 

If we write the functions G(9) and H (9) in the form 



G{9) = - exp (mXi - X x 9 + | v{&) + \-q(0)) 
H{9) = exp (-X 1 9 + 1 - |r/(0)) , 



(2.23) 



the functional relations (|2.19p and (|2.20p reduce to (|2.22p with a particular choice of the functions 
f{9) and g(9): 

i — m/3 cosh 0— 2ttiXq i „— mj3 cosh 9— 2iri\ 

/(0)=ln- g . fl ^ 9 rr -In (2.24) 

J v ' 1 — e -m/3coshfl+27riA _ g-m/3 cosh 6»+27r«A ' v 7 
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/ 1 p — m/3 cosh 9+2ni\o \ (-\ p — m/3 cosh 8— 2niXo \ 

g(6) = In — — . >{ J-. (2.25) 

M g— m/3 cosh 9-\-2iri\ g— m/3 cosh 6*—27riA j 

The branches of logarithms in (I2,24|) and ()2.25p are fixed so that for real 9 their imaginary parts 
lie in the interval (— it; it). 

The formulae (|2"7lTjh (I2T8I) . (f!T2"T|) and (l233l) - (|2T25l) provide a solution to the Dirac equation 
on the 1-punctured cylinder with specified branching (one should go back and check that all 
formal manipulations we have made with functions G{9) and H(8) indeed can be done). It 
remains only to verify the normalization condition (|2.16p . 

Let us take, say, the expansion ()2.18|) and rewrite it using the Poisson summation formula: 

CO 

' -e 6 



-mx cosh 9—im(y+kf3) sinh 9—2irik\ 



1 



The asymptotics of for \z\ — > is determined by the term with k = 0. Since Ai > 0, the main 
contribution to the corresponding integral is due to large \6\. Straightforward calculation shows 
that for Izl — > 



Mx,y) ~ ~ e- mX ^ 2+l ^ / 2 r(i - \ x )(mz/2)*- x l 

Z7T 



where 



oo 



\ r I \ g — m P cosh 9—2iri\Q -y g— m/3 cosh 8— 2ni\ \ 

i/oo = lim i/(0) = / d6 In - a — .. - In ^ . (2.26) 

0-»+oo 2lT J \ 1 — e ~mf3 cosh 9+2m\ y _ e -m/3 cosh 9+2ni\ J v ' 

—00 ' ' 

Therefore, the solution we have constructed differs from the element of the canonical basis only 
by a constant factor which can be set to unity by appropriate choice of A. Summarizing all 
these results, we obtain 

Theorem 2.3 The element of canonical basis on the cylinder with one branchpoint is given by 
the following expressions: 



« e ~\ 1 9 n+m xcosh9 n+ imysinh8 n I e ^ < Q ^ 

mpcosht'n VI/ 

w(x,J/)=A V — — g-Ai^-mxcoshen^m^/smh^ e ^ > Q ^ 

^— mpcosh0 n V 1 / 



neZ-A 

where the functions v{0), n(8) are determined from 



00 



ini I 1 —m/3 cosh 0'—27riAo i —m/3 cosh 9 1 —2-ni\ 

u(9)= I — tanh(0' - 0) In ■ .. - In . 

' 2lT \ 1 — e -mf3 cosh 9' +2m\ ^ _ g-m/3 cosh 0'+27riA 
^/ ^ g— m/3 cosh 0'+27riAo^j g— m/3 cosh #'—27riAo^ 



/df) 
— sech(9'-6)ln- 
2-7T 



^ _ e -m/3 cosh 0'+2ivi\ \M _ g-m/3cosh6»'-27riA \ 

— CXJ 

and normalization constant A = — 2sin7rAi e ~ w °°/ 2 . 
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3 Green function for the Dirac operator 

3.1 Green function for n = 

Let us calculate the Green function on the cylinder without branchpoints. The domain of the 
Dirac operator in this case consists of quasiperiodic functions ip(x,y + 0) = e 2mX °ip(x,y) that 
are square integrable in the strip S = {(x,y) : < y < (3}. After Fourier transformation 

oo 

= J^y J dxjdy i>(x,y) e~ l ^ + ^\ £, eR,^e y(Z + A ) 

-oo 

the Dirac operator and its inverse are represented by matrices 

D = -( m ~^ \ £>-! - 2 ( m %l 



-i£ m J ' m? + |£| 2 \ i£ m J ' 

where £ = £ x + i£<,, £ = £ x — £ K . Since the inverse transformation is given by 



f« -oo 

one obtains the following formula for the Green function Gq(x — x' , y — y') : 

oo 

U ^ (31) 

«f — oo 

Another two representations of the Green function will be useful for us. Choose, for example, 
x > and calculate the integrals in (|3.ip : 



p—mx cosh 9 n —imy sinh 6> n / i n \ 



/3 cosh l 

nGZ— Aq 



'C-UC+ 

Analogously, for x < one obtains 



d6 €~ mX cos ' 1 ®~ ^ m V sinh / j —(• ■" 
m /r- 2~7T 1 - e -«m/3sinh(9-27riAo ( -e"" 1 J ' ( 3 ' 3 ) 



3 mi cosh Q n +imy sinh 0„ 



<*><*.»> - e c ( x v j = < 3 - 4 > 

nGZ+Ao /" » \ / 

/" <i$ gWia; cosh 9+imy sinh ® / \ ffi \ 

= ~ m Jc-uc+ 2^1- e im P si * h9 - 2 ™ x ° \ e- 1 J ' ^ 3 ' 5 ^ 

In what follows, we shall also need the asymptotics of these expressions as y — > 0. To find 
it, let us rewrite (13.11) using the Poisson formula: 



oo oo 

Got-. y) = ^ E J J ( i{ I) s^r^p • (3-6) 
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The leading term of asymptotics is determined by the integral corresponding to k = 0, repre- 
senting the Green function on the plane. A little calculation shows that as \z\ — * 

m ( \n\z\ 1/z 

It is very convenient to write Green function as 

G {z) = 2G(z)J, J = 
since the rows of G(z) satisfy the Dirac equation (and not its adjoint, as the rows of Gq{z) do). 



i 
-i 



3.2 General properties of Green function 

The domain T> a ' X of the Dirac operator D a,x is chosen to consist of functions ip that have 
monodromies e 2wlXv (v = 0, ... ,n) and are integrable at |x| — ► oo. We also require (see (|2.6p ) 



{Z ^^^H^'H, v = l,...,n (3.7) 

where H 1 [a u ] denotes the space of functions that are square integrable in the neighborhood of 
{a u } together with their first derivatives. In the previous section we have shown that the Dirac 
equation D a,x tp = has no solutions in V a ' X . Thus "naively" we can think of D a,x as being an 
invertible operator. The kernel of the inverse is called the Green function G a,x . More precisely, 
the solution of 

D a ' X ifj = <p, 

is assumed to have the form 

ip(z) = J G a > x (z,z')Jcp(z')idz' Adz', (3.8) 

C\b 

Then one can try to determine the Green function by the following requirements: 

• The columns of G a,x (z,z') must satisfy Dirac equation D z G a,X (z, z') = for all z G 
C\(feU -{V}); they are square integrable functions at \x\ — * oo that have monodromy e 27nAl/ 
(u = 0, . . . , n) and singular behaviour (13. 7p at each singularity. Therefore, 

G a f(z,z')[a u ) = £ [4 V ] (z')w k+K [a u ) + b$ (z')w k *_ Xv [a u ]} . (3.9) 

• The integral operator with the kernel Dz' X G a,x (z, z') has to "cut out" the values of the 
function <p{z). Therefore, the singular behaviour of G a,x (z,z') for z — > z' must coincide 
with that of Green function for the Dirac operator on the cylinder without branchpoints, 

G a ' x (z, z') - G{z, z') eC l {z^ z'). (3.10) 

Remark. Suppose the function defined by these conditions exists. Then it is unique, since the 
columns of the difference of two such functions are obviously in W a '\ 

Next we determine how G a,x (z, z') depends on the second argument. To do this, let us define 
a matrix F a,x (z,z r ) satisfying the following conditions: 
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• The rows of F a > x (z,z') satisfy Dirac equation D z ,G^ X {z, z') = for all z' e C\{bl>{z}); 
they are square integrable functions at \x\ — ► oo, have the inverse monodromy e~ 2mXu 
(u = 0, . . . , n) and corresponding singular behaviour (|3.7p at the branching points: 

F;; X (z, z')\a v ] = [4 V J (z)w k -xAa v ) + 0$ (z)w k * +K [a„] } . (3.11) 

• The singular behaviour of F a,x (z, z') for z' — > z coincides with the asymptotics of "unper- 
turbed" Green function, 

F a ' A (z, z') - G(z, z') e C 1 {z' -> z). (3.12) 

Skipping the proof of existence of G a ' A (z, z') (or F a ' X (z, z')) we shall now prove the following 
Theorem 3.1 G a > x (z,z') = F a < x (z,z'). 

■ First we note an auxiliary relation: if f(z) and g{z) are smooth functions on some open set 
U CC, then 

{Df-Jg-f- JDg} dz A dz = {d z (f 2 g 2 ) ~ d s (fm)} dz A dz = d(f l9l dz + f 2 g 2 dz). (3.13) 
Now choose two distinct points x,y ^ b and 

(f(z)=F«:\x,z), 

\g(z) = G a ,f( z ,y), 

with U being the complement to the union of the disks f |J D e (a v ) \ U D e {x) U D e (y) and two 

strips: S' £ = {(x, y) : < y < e} and S" = {(x,y) : (3 — e < y < /?}. Integrating (|3.13p over this 
set and using Stokes theorem, one obtains 

= ^2 f Ui9idz + f 2 g 2 dz) + <j> {figidz + f 2 g 2 dz) + <j> {figidz + f 2 g 2 dz) . 

v=l dD E (av) dD e (x)U8D e (y) dS' £ UdS' s ' 

The expression under the integrals is single- valued on C\a, so in the limit e — > the last integral 
cancels out. The integrals over dD e {a u ) also vanish. One can easily check this, substituting 
instead of G a ' ? - (z, y) and F?l (x, z) their local expansions ()3.9p and (|3.1ip . It remains to calculate 

lim<! f F^[ X (x,z)Glj(z,y)dz + j> F^ {x, z)G%j(z,y)dz 



I 

and, similarly, 



8D E (x) 3D e (x) 

i-2m5iiG°£ x (x,y) - 2m5 i2 G^ x {x,y)\ = G^' x {x,y), 



a, A/ \ 7 / 7-,a,A/ \ / ~<a,\{ \ j— I itiO.A 



lim<j j F"{ {x, z)Gl' d {z, y)dz + j Ffy {x, z)G$ d (z, y)dz j> = - Jig {x, y). 

dD e (y) 8D £ (y) 
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Finally we get G a ^(x,y) = F?f(x,y). ■ 

Remark. Since the columns and rows of G a '^(z, z') satisfy Dirac equation in z and z' respec- 
tively, so do also the derivatives of the Green function d a)1 G a,x (z, z') and d afi G a,:K (z, z'). These 
derivatives are not singular at z — ► z' . However, their local expansions are more singular than 
(| 3 . 9 j) or (|3.1ip ; from the relations 



d z wi = YWi-i, \d. 



z w i = Y w i+v 



[d s wi = f wi+i, {dzwf = f w i-i 

it follows that 

d ait Gj' X (z, z')[a u ] = -y V«l/2j( Z ) U, -l/2-A I , + \ z ) w k-\ v + IkJ (^k+x,} > ( 3 - 14 ) 



fe>0 



d- afl Gf{z,z')[a v ] = -f (3-15) 



fe>0 

Therefore, if we somehow determine the coefficients ajj^jiz), Py^j^ 2 ^ anc ^ ^ nc ^ ^ e som ti° n OI 
Dirac equation with the same "extra" singular behaviour, it will coincide with the corresponding 
derivative of Green function. 

Let us consider a multivalued solution f(z) to the Dirac equation in the strip 5 = {(x,y) : 
< y < /?}, which is square integrable at \x\ — ► oo and has the following local expansions about 
the branchpoints a u (there are no other singularities!) 



In addition, only a finite number of coefficients CLu \ with k < is allowed to have non- 
zero values. Lower and upper continuations of f(z) across the branchcut d u (y = 0, . . . , n) are 

v . 

supposed to differ by monodromy multiplier exp(27ri ^ Afc) , so that the product Gfc (z, z')f(z'), 

k=o 

as the function of z', is single-valued on C\a. Using ()2.9p and Stokes theorem, one obtains 

™ f ^a,x f „ j 
2 J 

C\\jD s (fl„) 



G^.(z, z')f(z')idz' A dz' = (3.16) 



= E / G^{z,z')f x {z')idz'- j> G a ^{z, Z ')h(z')idz' = (3.17) 

" dD e {a v ) dD £ (z) 

= E / G^(z,z')M?)idz'+ j G a £{z,z')M7)idz'. (3.18) 

V dD E (a v ) 8D s (z) 

The surface integral (|3.16p does not converge for e — > 0. However, one can still compare the 
asymptotics of (13.170 and (|3.18j) for e —> 0: 

lim(l3l7D = -l^^(-l)fc+V2 sinwA ^H (2) ^_^. 2 ^) 5 
^ fc>0 
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lim (I335D = 1 E(" 1 ) fc+1/2 sin ^ a £j(*)b { -l + iSnM*). 

Finally we get 

EE(" 1 )' £+1/2sin ^ {P&H*)& k + a&H*j$\ = J^fa-m + SjiW)). (3.19) 

Now we shall make a special choice of f(z') to find the lowest coefficients of the Green function 
expansions. Analogously to the previous section (see (|2.12p ). let us introduce n special multi- 
valued solutions to Dirac equation, that are integrable at \x\ —* oo and have local expansions 

w M (A)[a„] = 5^ u w_ 1/2+ x v [a v ] + ^ {aV (™»( X )) w k+\»l a v] + b{ k\™»( X )) w k*-\A a v}} (3.20) 

fc>0 

The existence and uniqueness of such solutions follow from the existence of canonical basis; two 
bases coincide if all X u > 0. 

After the substitution f(z) = w M (— A) in (|3.19p one obtains 



^1/2,1 (*) ^ _ im ( W ^2(Z,-X) 



im 



■w*(z,-X). (3.21) 



V ^lAaW y 4 sm7rA M V w,a(z,-A) y 4sin^A M " MV 
On the other hand, the substitution f(z) = w* (z, A) in (|3.19p gives a formula for .{z), 



,1' = r wiz,A). (3.22) 

v / 4sm ^ 

Taking into account the earlier remarks (the formulae (|3.14p and (|3.15p ). we obtain an expression 
for the derivatives of the Green function in terms of the solutions (|3,20p : 



d aj G a ' x (z,z') = - m w j (z,X)^^ j (z',-X), (3.23) 

o Sill 7T/\jf 

d aj G a >\z,z') = -^l^w*(z,-X)®w*(z',X). (3.24) 

O Sin TTAj 

3.3 One-point Green function 

The formulae (|3.23p and (|3.24p can be used to calculate the Green function G a,x on the cylinder 
with one branching point {a}. Notice that in the case of a single puncture the solutions (|3.2U|) 
are easily expressed via the element of canonical basis. If we suppose for definiteness that 
< Ai < i then w(z,A) = w(z, A). Differentiating the local expansions, one can also verify 
that 

2 

w(z, —A) = — d z w(z, 1 — A), 
m 

Let us determine G a,x (z, z') if both z, z' G C\b are in the left half-strip: Re z, Re z' < Re a. Using 
the results of the previous section, one can write 

w(*,A) = -A(A) £ ^ j, (3.25) 
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*<„-A) = -A<l-A) £ G(8 -' A>e m , cosh()n ( ■ (3-26) 

Note the simple relations 

A[X)A(1-X) = 4sin 2 vrAi, 

v(0 ; A) = -u(9; 1 - A) = -u(9] -A), (3.27) 

r/(0; A) =t?(0; 1 - A) = 77(0; -A). 

If we substitute (|3.25p and (|3.26p into the equations f)3.23j) . (|3.24p and then integrate them, we 
obtain the following formula for Green function G a,x (z, z'): 

„a\, a x G(9i\X)G(6 n ;l - X)e e " 

G (z, z) — zsin7rAi £ £ ffl X (3.28) 

e f {( Z -a)e e i+(S-a) e - e i+( Z '-a)e e ™+(f'-a) e -«"} / e ^+e n ^ 



1 )+C(z,z>). 



The function C(z, z') does not depend on the position of branching point and has to be deter- 
mined. We remark that the double sum in the RHS of the last relation converges even if z = z r . 
This gives a hint that C(z, z') should be equal to "unperturbed" Green function G{z — z'; Ao). 
Assuming this is indeed the case, let us fix a = and rewrite (|3.28[) via contour integrals 



G 0,x (z, z') = imsmnXi 



de f do 1 G{e-\)G{e'-\-\)e e ' 



2yj- J 2tt (1 fAmfi sinh 9— 2-iriXo gim/3 sinh 9'+2ni\o ^ 

C-uc+ C-uc+ 



f {ze e +ze- 6 +z'e 6 ' +z'e- 9 ' } / 0+0' \ 

x IT^e' { e 0' 1 )+G(z-z';X ). (3.29) 

Analogously, if Rez,Rez' > 0, the Green function G°' X (z,z') is assumed to have the following 
form: 

d9 f dO' H{6; X)H(9'; 1 - X)e 9 ' 



G°' X (z,z')=imsm7TX 1 J ^ J 



2lT (\ — g— imf3sinh8— 2ni\\{-^ _ g— imj3 sinh 8'+2iri\\ 

C-UC+ C-UC+ 
„ -^{ze 6 +ze- g +z'e e ' +z'e~ e '} / 8+8' J 



x " ^8' { _ e e> 1 ) + G ( Z ~ *'■> A )- ( 3 - 3 °) 

In order to prove that the formulae (|3.29p and (|3.30p indeed represent the Green function, 
we shall construct their continuations to arbitrary values z, z' G C\b, and show that these 
continuations coincide with each other. 

Let us start, say, from the representation (|3.29p . At the first stage, we construct its contin- 
uation to arbitrary values of z only. This can be done by shifting the contours C_ and C+ in 
the integral over 9 to Im 9 = — ^ , Im 9 = ^ respectively. After this procedure one obtains 



G°> x (z,z')-G(z-z';X ) = imsinvrAi f — [ — 

J 2n J 2i 



d9 I d9 f G(9 f ' 1 \^^9'+mx' cosh. 6' +imy' sinh6' 



~^ gim/3 sinh 8'+2ni\o 

oo C-UC+ 
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-ie e+e ' 




-ie 


e 8 ' 




1 


ie e+e ' 


ie e 


) 


e 9 ' 


1 





(~}{Q Z7r/2' A)e — ^ mx s i nn S+my cosh 6 

' (1 _ e m,3 cosh 0-27riA o )(_f e ZjZ gfl') 

G(6> + itt/2; A) e im;rsinh9 - m ^ coshe 
^2 — e —mpcosh8—2mXo\U e 6 _|_ e 0') 

We can not do the same thing the second time, since the function standing in the integral over 
8' has the poles 8' = 8 ± However, one can shift C_ and C+ to the contours C e _f and C £ jf 
shown in the Fig. 4, and then let e — ► 0. The continuation of (13. 29ft to all z, z' 6 C\b is then 

oo . 

d# f d9> G(8+ l -f;X)G(8' + ^;1-A)e e 

■ X 



- A . ,. . . x ^ f d9 f dff G(8 + ^;X)G(8' + ^;1-X)e° 
G ' ( z ) = ?ra sm 7rAi > / — / — -. 3 — , „ n . . . . 5 — , „, , = ., . 

v ' ' 1 / / y 27r _/ 27r (1 — g— cm/3 cosh 0—27riAo)(l — e — fT " 1 ' 3 cosh e +2ni\ ^ 

+ (3.31) 



p am{ix sinhO—y cosh 0+ix' sinh O'—y' cosh 0') / „9+0' —A„8 



oo oo 

-4jnsiH7T/ 



oo . 

t» „ /" d0' G(0 + ^;X)G(8' - if-; 1 - A)e 9 



/-t, J 2ir 7 2ir 1 



^2 g— crm/3 cosh 0— 2ttiAo V J gtrm/3cosh 0'+27riAo ^ 

p am{ix sinh 0— j/ cosh 0— ia;' sinh O'+y' cosh 6') / „0+9' rr^cfi \ 1 



CT=±1 -oo -oo 



X ■ 

-cne" 



■ \ i 



The last two terms are "pole contributions" that can be calculated using (|3.27p and contour 
representations of the Green function on the cylinder without branchpoints. 



/ 

7T/2 


\Im9' 


2e 




c + 


w 







-^e > 








-71/2 


— >- 




Fi 


s. 4. 





Re 9 



If we perform the same manipulations with the representation (|3.30p in the right half-strip, 
the final expression will coincide with (|3.3ip due to the relations (|2.19p and (|2.20p satisfied by 
the functions G(9) and H{8). Thus the formulae (|3.29p ~ (|3.3ip indeed define the Green function 
G°' x {z,z'). 

4 Tau functions 

In this section we study the spaces of boundary values of some local solutions to the Dirac 
equation. These spaces can be embedded into an infinite-dimensional grassmannian. The r- 
functions are determined via the trivialization of the det*-bundle over this grassmannian. 
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4.1 Subspaces W int (a) and W ext (a) 

1/2 

Let us consider a circle L XQ = {(x,y)£C : x = xq}, and denote by (L xo ) the space of 

C 2 -valued quasiperiodic functions on L xo . Namely, if g G H 1 ^ 2 (L xo ), then g(y + /?) = e 2mX g(y). 
After Fourier transform the function g can be written as 

9(v) = Jf E ^n)e imi/sinh9 ", sinh0 n = -^n. 

1/2 

Let us introduce two operators, Q+ and Q_, acting on flY (L xo ) in the following way 



2vr 



Q±g{v) = j E Q±(^)5(^n)e^ sinhe " 



These operators have the properties of projectors, 

Q+ + Q_ = l, <£ = Q +J Q 2 _ = Q-, 

and thus define the splitting hV 2 {L X0 ) = © with iif^ 1 = Q±H 1 ^ 2 (L X0 ). One may easily 
verify that 



E IIQ±(W(^)|| 2 co S h0 n = | £ |5±(^)| 2 , 

j 

where 



2 

igZ+A ngZ+A 



$-(*») y ~ V ~Vn 1/2 v]l 2 J V 92 (On) 

Therefore, the function g is expressed through its polarization components g±(8 n ) as 

Let us now show that the elements of H7 (Ht) represent the boundary values of quasiperi- 
odic solutions to the Dirac equation in the right (left) half-strip x > xq (x < xq). To do this, 
rewrite the Dirac equation in the form 

d x rl> = ( ~ l&y ™ ) V- (4.2) 

\ m lOy ) 

If we put the initial condition ip(xo,y) = g(xo,y) with g(xo,y) E (i.e. all g+(0 n ) = 0), the 
solution of (14. 2j) in the right half-strip is 

_ 27T ^ e imy sinhe n -m(x-x ) cosh 6 n / ^1/2 _ w -l/2 \ / fj \ 

The convergence of this series is guaranteed by its convergence for x = xq. The solution of (|4.2p 
in the left half-strip can be constructed from the element of Ht in a similar fashion. 

The Green function G(z — z';X) on the cylinder without branchpoints provides a useful 
formula for the projections Q±. 
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Proposition 4.1 Consider the map Q : H^^L^) — ► H X {C\L XQ ) , defined by 



(Qg)(z)=i / G(z-z';X)a z g(y')dy' 



1 




[ ) , 5 e F 1 / 2 ^). (4.3) 



T/ien i/ie boundary values on L XQ of the restrictions of (Qg)(z) to the left and right half-strip 
are equal to Q+g and —Q-g respectively. 

■ To prove the proposition, one has only to substitute in (j4.3j) the Fourier expansion (|4,ip of g 
and the representations (|3.2|) and (|3.4|) of the Green function. ■ 

Suppose no two branchpoints have the same first coordinate. Then one can isolate the 

n 

branchcuts b\, . . . ,b n in the open strips Si,...,S n (Fig. 5). The union \J Sj will be denoted 

i=l 

by S. Consider the subspace of i? 1 (C\5) consisting of functions if) that satisfy on C\S the Dirac 
equation and appropriate quasiperiodicity conditions: 



tp(x,y + P) 



tp(x, y) for x < Xi, 

k 

exp{27ri J2 ^j}i , i x -,y) f° r xf<x< x% +1 , 

3=0 
n 

exp{27ri ^-^(a;, y) for x > x^- 
j=o 



We will denote by W ex t the space of boundary values of such functions. It is a subspace of W, 
the space of all quasiperiodic C 2 -valued H 1 / 2 functions on dS: 



r l/2 



1/2 



(4.4) 



Analogously, W{ n t C is defined as the space of boundary values of functions g € T> a ^ that 
solve D a ' x g = on S. 



!a R s, ! h, 



A 



• • • \y 



A 



x 

©> 



Fig. 5. 



The construction of the infinite-dimensional grassmannian in the next subsection heavily 
relies on the transversality of the subspaces W ex t and Wi n t in W . We postpone the proof of this 
fact; instead, let us explain how one can find the explicit formulae for the projections on these 
subspaces. Consider the restriction = g\* L o., ,aRq. of an element g E W to the boundary of 
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the strip Si. It is convenient to introduce the notation 

'-(*)•(*)■ 

(i) 

Example. Assume for a moment that the strip Si contains no branchcuts at all, i. e. both g L 

(i) 

and g R obey the same quasiperiodicity conditions, say, 



Then the map 



gf{y + /3) = e 2 ^gf{y), g® (y + (3) = e 2mX g$(y). 



Qg®(z) = i J G(z-z';X)a z g^(y')dy' = 
d L SiUd R Si 

G. t i(z - z; A)s? ) (z')dz' + G. )2 (z - z'; \)gf (z')dz' (4.5) 



d L SiUd R Si 



defines a function that satisfies the Dirac equation in Si- After a simple calculation, involving 
Fourier representations of g^' and of the Green function, one obtains the explicit formula 



, ,\ 2lT p m(x-x R ) cosh 9 n +imy sink 6 n ( ,1/2 _„- 1 / 2 \ / „(*) (a \ 

2vr „ e -m(x-xf)coshe n +imysinhe n / ,1/2 ^"V 2 \ / q \ 



+ 



Passing to boundary values, we see that Q induces a map on W. It is given by 



Q 



„(*) \ / _(*) \ / _(*) \ / n - \ / „(*) 



where {Cjg){6 n ) = e ( « « ^ g(0 n ) in Fourier representation. Furthermore, Q is a projec- 
ts) (i) 

tion onto the space of solutions to the Dirac equation on Si. If g L and g R represent boundary 
values of a function / that belongs to this space, the one-form in (|4.5p is closed, so the contour 
of integration can be shrunk to a small circle around z. Using the asymptotics of the Green 
function at z' — > z, one obtains Qf(z) = f{z). 

The generalization of the example we have just considered to the strip containing a branchcut 
leads to the main technical result of this subsection: 

Theorem 4.2 Suppose that G a,x (z,z') is the one-point Green function^ found in the previous 
chapter. Suppose that a £ S' , S' = {(x, y) £ C : xl < x < xr}. Consider the function g on dS' , 
which satisfies g\ gLS / £ H 1 ^ 2 (d L S') , g\ d n s , G H~^ 2 (d R S'). Then the map 

P s ,{a)g{z)= f G a $(z,z') gi {z')dz' + G a ${z,z')g 2 {z')dz> (4.6) 
d L S'Ud R S' 



2 Here a denotes a single point, and not the collection (oi, . . . ,a n ). I hope this abuse of notation will not 
confuse the reader. 
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defines a projection onto the space of functions f € T> a,x that solve D a,x f = on S' . The 
induced map of boundary values is determined by the following formula: 

™-(ZHZh(ZH$ !)(£)■ (47> 

where 

Onj r7r \ AAi + i -m(xij-a a; )(cosh6' ; +coshe n )-imay(sinh6';-sinh6» n ) 

(dg)(0j) = g 1 V — x (4.8) 

neZ+A 



2g-7riAi g j n _\ p—m(x R -a x ) cosh9 l +m{x L -a x ) cosh9 n —irna y (sinh9i—suih9 n ) 

= ^ £ ^ x («) 



_Ai + ^ — Ai + | 



vi-v n 



2g-?i"«Ai _< \ e m(x£— an:) cosh^j— rn(x_H— an:) cosh^n— imo y (sinhSi— sinhS„) 

W9)m) _ — E («») 

neZ+A 



x « t 1+! W +2 e ^(^-^)+|(^-^) (gw) / G z + Ao, 



- Ai + i m(xL ~a x ) (cosh 0; +cosh 9 n )—ima y (sinh — sinh 9 n ) 



and vi = v{0i; A), rji = 77 (0j; A). 

■ The derivation of (|4.2p - (|4.1ip can be carried out analogously to the previous example, using 
two more (in addition to (|3.28lJ3T3"0l) ) representations of the one-point Green function: 



n a,\, r\ ■ ■ \ sr^ G(9i; X)H(9 n ;l - X) e 
G ' (z,z ) = ^sln7^Al > > -k — — 

v ' ^ m R 2 rnsVi ft, rnsh ft.. 



9n+m(x—a x ) cosh 9i+im(y—a y ) sinhS; 

X 



' _ mj3 2 cosh 9[ cosh 6 n e° l — e dn 

ZeZ+Ao „ e z+A 



xe -m(x'-a x )cosh8 n -im(y'-a y )smh8 n I C « e« x fa. x < „ < J 



G^(z,z')=i S mnX 1 V V G(-6i,X)H(-d n ;l-X) e_ 



(4.12) 

e 0n 1 

—9 n —m(x—a x ) cosh 9i+im(y—a y ) sinh 0; 



m/3 2 cosh cosh # n e e ' — e 9n 

(4-13) 

^ e m(x'— a x ) cosh n —im(y'—a y ) sinh n / — e ' " — e ' j £ Qr x > a ■ > x' 
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When one applies the Stokes theorem to prove the projection property, the contour of integration 
in (|4,6p can be deformed into two small circles, around z and a. Using the expansions (|2.7p . 
(|3.1ip of the multivalued local solution to Dirac equation and Green function, one easily shows 
that the second integral vanishes. ■ 

1 /2 

Remark. Let us choose in HJ (L) a complete orthonormal family {ifk}, say, 

With these functions, we can find Schmidt norms of a, (3, 7, 5 and show that they are finite. 
For example, 

= K/Wo,¥V+A>| 2 = E E »A)| 2 , 



where f3(0i,6 n ) denotes the "kernel" of (3. This sum rapidly converges due to the exponential 
factors e -m{*R-a»)cabOi and e m(x L -a x ) COB he n in 0{e h e n ). Note, however, that in the limit p -» 
00, when a, (3, 7, 8 become integral operators, and 7 no longer belong to the Schmidt class 
due to the singularities in the kernels. 

We briefly outline the proof of the transversality of the subspaces Wi n t and W ex t in W, 
closely following Palmer's work [ID]. Suppose that / G W is decomposed as / = g + h, with 
g G Wi n t and /i G Wext- Theorem 12.11 guarantees the uniqueness of this decomposition, since the 
elements of Wi n t H Wext represent the boundary values of the functions from W a '\ It remains 
to prove only the existence. 

In order for g to be in Wi n t, one should satisfy the conditions 

(4.14) 





where dj, fa, ji, 5i are obtained from a, (3, 7, S by the substitution 

i— 1 i 
XL^xf, x R ^xf, a^ai, A ^E Afe ' A^^Afc. 

fc=0 fc=o 

Another set of relations follows from the assumption that h G W ex t- Indeed, should represent 
the boundary value of a solution to Dirac equation in the left half-strip x < x\ ; analogously, 
is the boundary value of a solution in the half-strip x > x^. This leads to two relations, 

4 1 ,- = 0, h$ + = 0. (4.15) 

(i) 

Next, the example we have considered shows that the functions h R and h L are boundary 
values of a solution in xf < x < xf +1 , if 

=u i h i i ) _, h ( § + =u i hf^\ i = l,-..,n-l, (4.16) 



where &i is obtained from uj by the substitution 



A:=0 
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One can transform (|4.15j) and (|4.16p into the conditions on g. Using (|4.14p to eliminate all g R 



(0 



(i) 

and g L + , we obtain a system of equations, 

5i _ - ^5^,+ + Ml-) - Jl- UiJR-i 



i = 1, . . . , n — 1, 
i = 1, . . . , n — 1, 



(i) _ f w w _ f 



(n) 



F (n) 



(4.17) 

(4.18) 
(4.19) 



If we introduce the notation 

Ui = - 







Vi 







i = 1, . . . , n — 1, 



#7 



Ak) r Ak-i) 



/ f (l) . ,(2) \ 

the system ()4. 17[) — (|4. 19[) can be rewritten as 



/ 1 *7i 
Vi 1 tf 2 
V 2 1 



;' = l,...,n, fc = 2, ...,n-l, 



An) 
JR,+ 



An) r An-1) 
J L - ~ Wn-lJn- 



\ f h\ ( Fi \ 




\ 



• £4-i 
V n -i 1 / 



V 9n J 



(4.20) 



The operator, standing in the LHS of (|4.20p . represents a compact perturbation of the identity, 
and thus is Fredholm of index zero. Since to every nontrivial element of its kernel there corre- 
sponds a nontrivial element of W{ n t D W ex t, the kernel should be zero. Therefore, this operator 
is invertible and can be used to construct the decomposition / = g + h explicitly. 



4.2 Grassmannian, det*-bundle and its trivialization 

First we introduce several important definitions, following Segal and Wilson [18] (further details 
can be found in [T4"l [19]). 

Definition 4.3 Suppose we have a complex Hilbert space H with a given decomposition H = 
H + © H-. The Grassmannian Gr(H) is a set of all closed subspaces V C H such that 

• the projection pr + : V — > H + along H_ is a Fredholm operator; 

• the projection pr_ : V — > H- along H + is a Hilbert- Schmidt operator. 

The first requirement means that the codimension of V H H+ is finite in both V and H + . The 
connected components of the Grassmannian are distinguished by the value of the index of pr + . 
We shall work only with the component Gro(H) that corresponds to zero index. 
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Definition 4.4 The invertible linear map v : H + — > V is called an admissible frame for the 
subspace V £ Gro(H) if pr + o v : H + — > H + is a trace class perturbation of the identity. The 
fiber of the del* -bundle over V consists of the equivalence classes of pairs (v,a), where v is an 
admissible frame, a is a complex number and (vi,ai) ~ {1)2,(22) if ot\ = 02 det . The 

canonical section of the det* -bundle is given by a : V 1— ► (v, det (pr + o v)) . 

In the work of Segal and Wilson, as well as in almost all subsequent papers on the subject, 
the Hilbert space H is the space of all square-integrable complex-valued functions on the unit 
circle S 1 = {z £ C : \z\ = 1}. H + and H- are spanned by the elements {z k } with k > and 
k < 0, respectively. 

We are interested in a more complicated model of the Grassmannian. H is identified with 
the space W (see (I4.4h ) of square-integrable quasiperiodic C 2 - valued functions on the boundary 
dS. Let us fix a collection of points a = (a°, • • • , a°) such that a® G Sj, j = 1, . . . , n. Then one 
can define the Grassmannian Gr{W) with respect to the splitting W = Wi n t(a°) © W ex t- The 
crucial observation, similar to one made by Palmer in [TO], is that Wi n t(a) £ Gtq{W). 

Now introduce two admissible frames for the subspace Wj n t(a). The first one, which will be 
denoted as P(a) : Wi n t(a°) — > Wj n t(a), represents the projection of Wi n t(a°) on Wi n t{a) along 
W ex t- It is easy to guess that P(a) inverts pr + . Thus the canonical section can be written 
as a : Wi n t(a) 1— » (P(a),l). The second admissible frame, F(a) : Wj n t(a°) — > Wi n t(a), is the 
restriction to Wi n t(a ) of the direct sum of the appropriate one-point projections (|4.6p ~ (|4.7p : 

F(a) = (P5 1 (ai) © . . . © Pg n (a n )) . It defines a second, trivializing section 1? : Wi n t(a) 1— > 

Win* (a ) 

(P(a), 1). The determinant of the Dirac operator D a ' A , or the r-function, is determined from 
the comparison of the two sections, 

r(a, a°) = = det (P(a)- 1 P(a)) . (4.21) 

v(Wi nt {a)) 

Remark. This ideology originates from the work [S], where the isomonodromic r-function, 
associated to a fuchsian system on CP 1 , was interpreted as the determinant of a Cauchy-Riemann 
operator, whose domain incorporates functions with specified branching. 

In order to calculate r(a, a ) more explicitly, let us use several results and notations from 
the previous subsection. Suppose that / G Wi n t(a°), g £ Wi n t(a), then 

r j) i.^A'nir u j: //r-\> ; ,, ; . w = (^ I'^j)- ( 4 - 22 ) 

The functions / and 5 can be represented by columns 

/ = (fl ■ ■ ■ fnf, 9 = (ffl • • -gnf- 

The map P(a) : Wi n t(a°) — > Wi n t(a) is given in these coordinates by the identity transformation. 
To obtain the representation of P(a) _1 , for each g £ Wi n t(a) one should find a function / £ 
Wi n t(a ) such that g = f — h with h £ W ex t- This amounts to the same calculation as we 
have done earlier in (|4.15p - (|4,20p . Taking into account the additional condition (|4.22p on /, one 
finally obtains 

(1 + M(a))g= (1+M(a ))/, 
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where 



M(a) 



( Ui{a) 
Vi{a) U 2 (a) 
V 2 (a) 



V 

Therefore, the r-function is equal to 



\ 



U n -i(a) 
V n -i(a) / 



r(a, a ) =det{(l + M(a)) (l + M(a )) X } . (4.23) 
In fact, one can derive an even more convenient representation. Let us introduce the matrix 



M{a) 



( Ui(a) 
Fi (a) C/ 2 ( a ) 
f 2 (a) 



V 



\ 


U n -i{a) 
V^-i(o) 



with 



17,- 



-«2>i7j+i 






-tyfy 



j = l,...,n- 1. 



The matrix 1 + M(a) is a product of an upper triangular and a lower triangular matrix with 
identities on the diagonals. Thus we have 



det<^ (l + M(a )) (l + M(a 



-i 



1. 



Multiplying the RHS of (|4.23[) by this determinant, one finds 



-(a, a ) = det { (l + M(a))'\l + M(a)) (l + M(a )) 1 (l + M(a )) } = 



where 



r a 



) = det { (1 + M(a)) 1 (l + M(a))}. 



(4.24) 



Example. Consider the simplest nontrivial situation, when there are only two branching points 
on the cylinder. Let us introduce the notation 

A = Ao + Ai, A = Ao + Ai + A2, a x = (a,2) x — {a-i)x, a y = {a2) y — {a\)y 

In the case n = 2 the inverse matrix (l + M(a)) 1 looks particularly simple, 

(1 + M(a)) _1 = 



1 -Ux(a) 
-Vx{a) 1 

Using this formula, one can find that the two-point r-function is given by 

r(a) = det (1 - K) , 



(4.25) 
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where the operator K = Q}\ct\Cji52 can be represented by the infinite-dimensional matrix with 
the entries 

K = 4 sin 7T Ai sin 7rA 2 (v m v n ) 2 x {^26) 
/3 2 y/ cosh 6 m cosh 6 n 

, . „ I cosh 6 m +2 cosh fy+cosh n . sinh 6 m — 2 sinh S ; +sinh S n p(6 m ) + 2p(6 l )+p(6 n ) 
^Ai A2-ri-g — TTlla^l 2 h^^Cty 2 1 2 

X x 



(1 + v m uj)(l + viv n ) cosh (9/ 
The indices take on the values m, n E Z + A and 

2p(0) = r/(0; A, A) - 7/(0; A , A) + ii/(0; A, A) - iu(0; A , A). 
One can also write K as 

K = 4sin7rAi sin7rA 2 • VV T , 

_ , > ^1 — ^2+1 I | cosh 9 m +cosh 9 ;1 , „ sinh 9 m -sinh 9 n , p(0m) + p(8n) 

T/ 1 fm^n 3 e -™l««l 2 2 + 2 - 

v mii = ^ , m,Ji£i; + A. 

p Vcosh0 m cosh6 , „(l + v m v n ) 

These explicit formulae for the r-function are in some sense a reward for the technical work put 
in the calculation of the element of canonical basis on the 1-punctured cylinder (Theorem 12 .3p . It 
would be interesting to compare them with the correlation functions of twisted fields, calculated 
in the lattice regularization of the Dirac theory on the plane [1] . 

Remark that the final answer (|4,25p ^ (|4.26p for the two-point r-function is independent of 
the choice of localization (coordinates of edges of the strips Si,... ,S n ). To show that this is 
true in the general case, we prove 



Proposition 4.5 The logarithmic derivatives of the r-function \^.21 ) are given by 



dlnT{a,a°) = ^J2) °i/a da » + a 1/2 CM -A)) da v \ . (4.27) 

■ Consider the n-point Green function G a,x (z, z') and construct the map P{a) : W — > Wi n t(a) 
in the following way: 

P(a)f(z) = J G a }(z,z')f x {z')dz' + G a ; 2 \z,z')f 2 (z')d?. (4.28) 
dS 

It is easy to see that this map defines the projection on Wi n t(a) along W ex t- Indeed, let us 
write the function f(z) as / = g + h, with g G Wi n t(a) and h S W ex t- The form in the integral 
P(a)g(z) is closed, thus each piece dS^ of the integration contour can be shrunk up to two small 
circles, around z' = z and z' = dp. Computing the residues, we obtain P{a)g{z) = g(z). In a 
similar fashion one also shows that P(a)h(z) = 0. 

It is clear that the admissible frame P(a) : Wi n t(a°) — ► Wi n t(a) and the projection pr + : 
Wi n t(a) — > Wi n t(a ) are the restrictions 



P(a) = P(a) pr + = P(a c 

W int {a°) 



W int (a) 



Let us analogously consider the map F(a) : W — ► Wi n t(a), which is by definition the direct 
sum of the one-point projections, F(a) = Ps 1 («i) © . . . © Ps n K)- The second admissible frame 
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. Its 

W int (aP) 



Wi„t{a) 



that we have used, F(a) : Wi n t(a°) — > Wi n t(a), is the restriction F(a) = F(a) 

inverse, which we denote as F(aP) : Wi n t(a) — > Wj n t(a ), is given by F(a°) = F(a°) 
Therefore, differentiating (14.211) . one obtains 

din r(a, a ) = -Tr {d (F(a) -1 P(a)) P(a) _1 F(a)} = -Tr {F(a°)d(P(a))pr + F(a)} . 

Recall that the traces in the last formula, and the determinant in (|4.2ip as well, are calculated 
on the subspace Wi n t(a°). However, since the range of both P(a°) and F(a°) is Wj n t(a°), we 
can forget about this restriction and replace under the last trace P(a) by P(a), F(a) by P(a), 
pr + by P(a ) and F(a ) by F(a ). If we also use the relations 

P(o)(l - P(a )) = 0, P(a)(l - P(a )) = 0, 

the trace becomes 

dlnr(a,a°) = -Tr w {f (a ) d{P (a)) P (a )F (a)) = -Tr w [d(P(a))F(a)} . 

Taking into account the explicit form (|4.28p of P(a) and the formulae (|3.23p - (|3,24p for 
the derivatives of the Green function, one can show that dP(a) is an integral operator with 
degenerate kernel. Then we have, for example, 

2 n f 

a a ,lnr(a,a°) = ^- x -^ J {(w„(s,-A)) ^(a^z, A)) dz+ 



dS u 



w„(z, —A) J 2 (p Stt yw„(z, A)j dz J . 



Again applying the Stokes theorem, each contour dS^ can be shrunk to a small circle around 
a«. Only one circle, around a„, gives a non-zero contribution, which can be calculated using the 
asymptotics of the one-point Green function. At the end of this calculation one finds 



d av lnr(a,a°) = y afy(w v (\)), d &v In r(a, a ) = — af^w^-A)), 
as claimed. ■ 



5 Deformation equations 



Let us now find the differential equations satisfied by the elements (I3.20p , using the idea that we 
have already exploited in the calculation of the derivative of the Green function. For example, 
consider the solution w^(A) and differentiate it with respect to a p . We obtain again a solution 
of the Dirac equation but with more singular local expansions at the branchpoints: 



k>0 



m 

"2 V 



l + EKW*)) Wk-t+M +6 ( r } (w M (A)) wt +l _ Xv [a u }} 

fe>0 
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Adding appropriate linear combination of {w^(A)}, {d z w v (X)} and {dfW^A)} (rj = 1, . . . , n), 
one can annihilate the coefficients near the "extra" terms W-s/2+\ v an d ^-i/2+a„- Then the 
result will vanish identically, since it is clearly in W a ' A . This observation can be written in the 
following general form: 

da,aW(A) = ($ 8 Z + <5>*8- z + *) w(A). (5.1) 

n 

Here d aB , = ^2 (dcij ■ d aj + dtlj ■ . ) denotes the differential with respect to the positions of the 
j=i 

singularities, <£, <3?* and ^ are matrix- valued one-forms, and w(A) = (wi(A) . . . w n (A)) T . 
Let us introduce the notation 



/x,v=l,...,n 



/x,i^=l,...,n 



In particular, one has Co = 1 and Cj = C? = for j < 0. The system (|5.ip . being rewritten in 
terms of {Cj}, {C*}, amounts to 

777 777 _ 777 777 

dCj - - C J+1 dA - - Cj-xdA = - $ Cj+i + - $*C J - 1 + (5.3) 

777 777 — 777 777 

dC/ - - CjUdA - - C* +1 dA = -$ Cj U + - **Cf +1 + 9Cf, (5.4) 

where dA = {S^daj,)^^...^ and dA = (^da„) Mli ,=i,...,n- 

Note that the expansion coefficients obey a set of algebraic relations. To derive them, let us 
first consider two multivalued solutions to Dirac equation, u and v, that are square integrable 
at | a: | — > oo and have the local expansions (|2.3p at the singularities. We shall assume that there 
exists a negative half-integer number kg such that 

M ( n ) = («) = («) = («) = 0, v = 1, . . . , n, 

for all k < ko. Using (|2.9|) and Stokes theorem, calculate in two different ways the integral 

j (u\V\ + U2V2) idz A dz = im <j> U2V\dz = —im <j> u\V2dz. 

C\\jD e (a v ) \JdD e (a„) \JdD e (a v ) 

Comparing the asymptotics of the corresponding boundary integrals as e —* 0, one obtains 

E E {^^IW-^^W^-l^sinvrA^O. (5.5) 

If we now put u = w M , tj = w p , then (|5.5p leads to the relation (analogous to (|2.13l) ) 



b i%(™f>) sin7r \> = & i/ 2 (w P ) sin 7^, 



or, in matrix notation, 



C Q * sinvrA = [C *sin^A] T , A = (^A„) ftl(=1 „ . (5.6) 

On the other hand, the substitution u = w M (A), v = w*(— A) gives 

Ci(A)sin7rA= [Ci(-A) sin7rA] T . (5.7) 



28 



Finally, observe that the entries of the n-dimensional vector <9 s w(A) — |C *(A)w*(-A) 
belong to W a,A and thus are all equal to zero. This gives two more relations, 



C *(A)d(-A) = C*(A), C *(A)C *(-A) = 1. (5.8) 

It is easy to see that for positive A^, \ u both w^(A),w l/ (A) G W a,A , so we can calculate the 
inner product: 



<w M (A),w„(A)) = -46f / ) 2 (w^(A)) sinvrA,. 

This shows that the submatrix of C *(A) sin 7rA, associated with the indices corresponding to 
positive {Ap}, is negative definite. In a similar fashion, one finds 



(w;(-A),wJ(-A)> =4&$(w„(-A)) sinvrA M 

for negative Au, A^. Combining this formula with the second relation in (|5,8p . one can prove 
that the submatrix of (C *(A)) _1 sin7rA, that corresponds to the "negative" indices, is positive 
definite. 

Let us return to deformation equations (|5.3p and (|5.4p . In order to determine the unknown 
matrix- valued forms $ and $*, let j = —1. One then finds 

$ = -dA, = -C *dA {Cjy 1 . (5.9) 

Specializing to the case j = 0, we calculate the form \E' and obtain a matrix equation, 

* = | [cL4, d] = dC * (C^)- 1 + f [C *dl (Co*)" 1 , d* (df)" 1 ] . (5.10) 

For j = 1, higher order coefficients arise. However, "antiholomorphic" part of (|5.3p and "holo- 
morphic" part of (|5.4p comprise only the coefficients that are already involved: 



Cq^CoT 1 =0, (5.11) 



d- a Ci + 

TJl 777 

d a C{ + - [dA, C *] - - [dA, d] d* = 0, (5.12) 

n n 

where d a = ^ daj • <9 aj and d s = X) ' ^% in addition, the diagonal part of (|5.3|) implies 
./ i ' i=i 

777 

4 diag Ci = y diag ([dA, d] d) . (5.13) 

In order to write the deformation equations in more compact and standard form, introduce 
the notation 

G = Co*sinvrA, = y [<LA, C x ], 9 f = 9 T . (5.14) 

Using the symmetry relations (|5.6p - (|5.8p . one can show that (|5.10p transforms into 

dG = OG + G@l (5.15) 
Instead of the equations (|5.1ip and (|5.12p we have two conjugate relations 

777 i _ 777 m m j 

d- a C 1 = -[dA,G}G-\ d a d = - [dA,G]G \ (5.16) 
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and the last equation (|5.13|) can be rewritten as 

d a diagCi = diag(eCi). (5.17) 

We easily find from (|5.14p and (|5.15p that det G = const. It is also very instructive to deduce 
the closedness of the 1-form 

n = ?f^{ a 1 /l(^(A))^ + 4/ ) 2 (w i ,(-A)) da v \ = ^Tr (C 1 dA + 'D 1 dA), 

standing in the RHS of (|4.27p . from the deformation equations. Indeed, 

dfl = j Tr (@ d A dA + j [dA, G] G' 1 A dA + G C 1 AdA + j [dA,G] G' 1 A dA 

2 

= -^-Tr [C\dA A C\dA + C~\dA A C\dA + GdA A G~ 1 dA + GdA A G ^dAj = 0, 
so the form f2 does represent the differential of a function. 

Example. As an illustration, let us find the explicit form of the deformation equations in the 
case n = 2. Suppose that Ai > and A2 < 0. Then G\\ < 0, det G < 0, and the matrix G can 
be parametrized in the following way: 

„ / — e^smip e llf cos ib \ , 

\ e lv cosip e v smip J nrir 

where < < ir, x > 0- We shall also denote 

C\ = ( ^ n ) , q = m(a 2 - oi)/2, g = m(a 2 - ai)/2. 



_ A21 A22 , 
From (|5.15p one obtains 

dG _ x ( Ai 2 e~ i</3 cosV' A^e^sm^ 



<9q X \ A 2 ie r 'sim/' -k2ie lip costp 

dG _ x ( Aue tip cosip A 2 ie v sin-0 
dq 

This leads to the relations 



dq X V Ai2e _,? sin-0 -A 2 ie~ ;i ' / ' cos V> 



(5.18) 



= i tg 2 ip d q r], dqip = -i tg 2 tp d q r], 

Ai 2 = e"+^ (5 g V - i ctg i/> d g ip) , A 2 i = e""-^ + i ctg V 3,?) • 
Next, the first formula in (|5.16p implies that 

dC\ ( cos 2 ^) e v+up cos ^ s i n -0 

dq \ e" 11 '^ cos ^ sin ijj — cos 2 ^ 

Off-diagonal part of this relation leads to a system of coupled differential equations, 

/ + iiP^ ^ ^ + sin V 7 cos V- = 0, 

l d lW = sin^cosV ^ ^ + ^ • 

Finally, the formula (|5.17p and the diagonal part of f|5. 18|) give the second logarithmic derivatives 
of the r-function: 

( dqqlllT = COS 2 tp, 

dqqlnr = (dqi;) 2 + ctg 2 i;(dq V ) 2 , (5.19) 

dqqlllT = (dq?p) 2 + CtgV (dqip) 2 . 
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6 Discussion 



When one tries to generalize the above theory, a natural question arises: if it is possible to 
develop the theory of monodromy preserving deformations for the massive Dirac operator on 
the arbitrary two-dimensional surface M with a metric? It appears that M should be then a 
homogeneous space for a group G, acting on M by isometries. There are only five such surfaces: 
plane (G = E(2)), cylinder and torus (G = T 2 ), Poincare disk (G = PSU(1, 1)) and the sphere 
(G = PSU(2)). The plane and hyperbolic disk were studied earlier by different authors (see 
references in the Introduction). The present paper is devoted to the cylindrical geometry. It is 
interesting to note that the derivation of all "implicit" results (factorized form of the derivatives 
of Green functions, deformation equations, etc.) can be transferred almost literally to the case 
of torus. What is even more important — the technical results, obtained in this work (namely, 
the formulae for the one-point projections in the Theorem 14,2ft allow to calculate the r-functions 
on the torus explicitly. I hope to discuss these matters in greater detail elsewhere. 

The second task is to give a proper formulation and solution of the problem in the quantum 
field theory language. Let us interpret the coordinate along the cylinder axis as time, with the 
space coordinate living on the circle. The time axis is splitted by the branchcuts &i, . . . , b n into 
n + 1 intervals. The evolution in each interval is governed by the Dirac hamiltonian, which is 
diagonalized in the free-fermion basis. These free fermions, however, obey different periodicity 
conditions (statistics) in different intervals. Corresponding Fock spaces are intertwined by the 
monodromy fields, whose correlation functions can be written in terms of Lehmann expansion 
over intermediate eigenstates of the hamiltonians. (In the two-point case, this corresponds to 
the expansion of the determinant (|4.25p ). The problem transforms then into the calculation of 
form factors of monodromy fields in the finite volume. 

Another important problem is the investigation of the ultraviolet (m — > 0) asymptotics of the 
r-functions. On the plane, the connection between the Ising model and singular Dirac operators 
was already used in [TTJ to give a rigorous proof of the Luther-Peschel formula. 
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